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ABSTRACT 



A mathematical model of a branching stochastic process utilizing 
generating functions is presented,, The probability distribution 
of the number of members of the process at discrete time periods, 
z n , and the probability of extinction is discussed,, When there is 
a non-zero probability of surviving indefinitely, the normed random 
variables, z n /E(z n ), converge with probability one; the cumulative 
distribution of this random variable is absolutely continuous. The 
time until extinction and the total number of members of the process 
is examined when the probability of extinction is one. The distri- 
bution of the z n given that z n is not zero is discussed for this case. 
The maximum likelihood estimates for the probabilities involved in 
the process are determined. An example is given of a branching pro- 
cess in which the probabilities are dependent on time and a solution 
is found using Laplace transform methods. 
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NOMENCLATURE 



Symbol 
P(A) 
P(A B) 
P(A,B) 

z r 

X 

A i 



f(8) 

Pr 

f n( s ) 

Pnr 



G n (u) 

n( s ) 

Q(u) 

k(s) 

H(u) 

(s) 



Meaning 

the probability of the event A 

the conditional probability of A given B 

the joint probability of A and B 

the number of members of the r^ generation 

the number of offspring of the member of the 

n^h generation 

the expected value of z^ 

the variance of z 

i 

generating function of z^ 

probability of r members in the first generation 
generating function of z n 

probability of r members in the n^h generation 
probability of extinction 

z 

a normed random variable ( 11 ) 

ri 

the limit of w n as n approaches infinity 
the cumulative probability function of w n 
the characteristic function of w n 
asymptotic distribution branching from f(s) 
probability generating function 

'V,' . 

the asymptotic distribution branching from k(s) 
the characteristic function of H(u) 
probability that z n ^ 0 

probability that z n equals r given that it is not 
equal to zero 



v 



y n number of members of the n^h generation given that it is 

not zero 

Fn(s) probability generating function of y n 
N number of generations to extinction 

(s) moment generating function of N 

S total number of individuals in the process 

g(s) probability generating function of S 

z mr number of individuals in the m^b generation who have 

exactly r descendents in the (m l) s ^ generation 
S n total number of individuals in the first n generations 
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INTRODUCTION 



The present paper is an essay on the theory and application of 
discrete branching stochastic processes,, Most of the asymptotic 
theory was extracted from the work of Harris [10] who derived many 
of the relationships for this class of stochastic processes,, Let 
us introduce these processes by means of a historical example,, 

In l87h, one of the first investigations into discrete branching 
processes was made by the Rev. H. W. Watson,, Dr„ F. Galton had 
posed the problem of extinction of family names to the mathematical 
community in England since many names, prominent in history, were 
no longer to be found. The problem was posed as follows: Assume 

a large nation with N adult males, each with a different surname. 

What proportion of names would be extinct after r generations? In 
the r^h generation, how many people would have the same name? Rev. 
Watson [8] derived some relationships that enabled him to compute 
the probabality of extinction and made some sample calculations to 
demonstrate the difficulty of computing numerical answers. Notable 
among these relationships was the functional relation of the proba- 
bility generating functions involved (section III, theorem 1). 

Branching processes received some attention thereafter and the 
theory was applied to problems of gene mutation [7], population 
growth (jL2j , cosmic radiation [l] , ionization fluctuations [l8l , 
neutron production in atomic piles [ll] , and cascade theory [20] „ 

In this paper we will review the theory of simple discrete 
branching processes and some of its applications. Let us consider 
the original problem of survival of family names and state the 
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basic problem of branching processes in that framework,, Suppose 
an immigrant by the name of Door comes to the United States and 
there is no one else here having that surname , What is the prob- 
ability that the name survives r generations and how many men will 
there be having that last name? 

There is a probability that Mr„ Door will either not marry or 
not have any male children; there are other probabilities associated 
with the number of sons he may have,. We will denote these probabili- 
ties by p r , i,e,, the probability of no sons will be p c , of r sons, 
p r . We need a symbol for the number of individuals with the name 
of Door in the n't* 1 generation, so let z n be the number of individuals 
with the name of Door in the nth generation,, Then we can write 
P(z 1 2 r) = p r indicating that the probability of r individuals with 
the name of Door in the first generation is p r , If Mr c Door coes 
have sons, each of them may have children and ti.e process is continued. 
This may be represented as a tree type graph with an initial node 
representing his sons, and so on„ A method of analysis of the 
problem based on graph theory was proposed by Otter [19] „ In this 
paper the method of generating functions is used. 

We will assume that the number of sons born to one of Mr, 

Door's sons follows the same probability law as that of Mr, Door, 

In other words, the probability that Door, jr„, has n children is 
p n „ Then, if is the number of offspring of the it* 1 individual 

in the nth generation, we have P(z n+ ]_ = rjz n s m) = P(^ ~ n) „ 

is J- 

This is the outline of the problem, the notation, and the assump- 
tions made in this paper. Before taking up the mathematical theory, 
some other problems will be posed in this framework. 
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In atomic piles, the number of free neutrons is an important 
factor. If we start with one neutron and a fission takes place, the 
probability that r neutrons are released is a nuclear constant which 
must be empirically determined for the amount and arrangement of active 
material. The probability of no fission (no descendents) is the 

average probability of leakage or absorption. An atomic pile in 

♦ 

a nuclear reactor is subcritical when the reactor is not operating 
and there is no possibility of an explosion, i.e„, the probability 
that the number of fissions increases sufficiently to cause an 
explosion is zero; the mean of the probability law for a subcritical 
mass is less than one. As the amount of fissionable material is 
increased, the probability of an explosion finally becomes greater 
than zero. When this occurs, the mass is termed supercritical and 
the mean of the probability law for this situation is found to be 
greater than one. A nuclear reactor is termed critical When it is , 
operating; the mean of its probability distribution is exactly one. 

The organization of the paper is as follows: The general math- 

ematical character of a discrete branching process is presented in 
section I. It will be seen that the mean number of progeny per 

individual plays an important rol,e in the asymptotic theory. This 

% 

theory is presented in section II for the mean of the probability 
distribution greater than one and in section III if it is 'not 
greater than one. The maximum likelihood estimation of the distri- 
bution of progeny per individual is given in section IV. An example 
of an ionization cascade is presented in section V. Tlds is a more 
general branching process than those treated in the model used in the 
earlier sections and its treatment suggests how the methods of Markov 
processes may be applied. 
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GENERAL 



The structure of branching processes was presented in the 
introduction where several assumptions were made and some notation 
was introduced. It is appropriate now to summarize the assumptions 
made in this paper and to extend the notation. 

Consider two related classes of random variables: 

z n the number of "individuals” in the n^' n generation 
the number of "sons" of the i th individual in the 
n th generation, 

They are related by the equations: 

Z h - Xl for n= 1,2,3,. . . 

The quantity z Q is the number of the initial group who start the 
process; z Q is assumed to be one, since this does not affect the 
generality of the theory as is pointed out at the end of this section. 
The probability distribution of the first generation is given by p r , 
i.e., the probability that there are r members of the first generation 
is p r . 

The following assumptions are made throughout this paper: 

(1) Given z n is equal to m, the random variable, z n ^q , is the sum 

of m independent random variables, each one of which has the same 
distribution as z^ : = f) = Pr 

(2) The variance of z ^ is finite. 

(3) P r / 1, for all r ; otherwise the process degenerates to a simple 
arithmetic calculation. 

(U) p 0 + P]_ / 1, otherwise we have the binomial case whose treatment 
in .the context of branching processes is trivial. 
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The following notation is used throughout: 



the mean of Z]_ 
the variance of z^ 
the probability generating function 

of Z]_ (s is a complex variable) 

p„ r = P(Zn = r) 

oo 

Pnr the probability generating function 

of z n . 

From the assumptions listed above, it is clear that : 
f Q (s) = s, p^ r = p r , f'(s) and f"(s) are continuous on the set of 
points consisting of |s|<l and s=l c 

Several authors (see [2], [9], jjll] ) have examined the generating 
functions and have found some basic relationships. Some of these are 
presented below. 

Theorem 1 . f n (s) -is the n^* 1 functional iterate of f(s) : 



A = E.Ui)-S rp 

/ r=o ; 

<r^Vor(? ( ) = S(r A p r \-^ 

Rs) = ]> p r s r 



F n (s) - f (f nH (sJ) - Fn-i ( F ( 



Proof : From the initial assumptions, the are independent for all 

1 and n > P (Xi'= k) = p k z n = 2* xT‘ ,J 

Ul 



It follows that PUn+i -K|z n =j) 



is the coefficient of s^ 



in the power series expansion of (f(s))^ since the generating function 
of the sum of j mutually independent random variables is the product 
of their generating functions. See Feller [6, p. 2!?0ff.] „ 
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PO I 

Hence, = 2 PjM = f(f(S« 
F 3 '5)=2p u [fls)] J = Uffsfl =FCf(f(sj)3 -fif^isfl 

^ J-0 J 



0O 



f n (s)-^ p nH [f(sj] J =- F n . a cf(s)l L . (f(f(s)J) . . .)) 

j * o 

T 

- f(f n J,(Sl) 

«• 

The coefficient of s r in f n (s) is the probability that there 
* 

are r individuals in the n^* 1 generation. Thus the questions posed 
by Dr. Galton may be answered in the framework of the present model. 
The constant term, p no , is the probability that the family name is 
extinct by the n^h generation. 

The mean and the variance of z n can be determined using theorem 1. 
Theory. EUn)=/U n U 



- rus 



! J - 






CO 



— F> = ^r&, r ~ E UnJ 



f n (i) = 2r(r-ij p nr = E (z n ) - £ lz„) 



|| | ^ 

Therefore, V<ar ( 2 p, ) — 0) + f^(l) “ [£>,'(*)] , 

By theorem 1 ( S) = Cf ($)] — f Cf M (S)J 

Taking derivatives, f^j(s) = ^ [f(s)J^ (S.) — f C f n (s)] (S) 
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Fn'.lSl = f" Cf (s)] [f'lsO^ t f'[f (S i]f“(s) 

- F"[f n (S)3Cf„'(Sfl i + f'Cf n WJK'M 

Hence (ll = /U £(2n| - M* 1 

Pn CFdllCF'HU 4 + f„' [f ^ 



Solving for fj! t (1) gives 

n ' [f'UIJ 1 - f'(l) 

/A.* - /U 



g a yU, n 



- yU. n tyU' 






Thus Var ( z n ) = £■ - - 

/Ur 



for Li 



For the case^sl , L'Hopital's Rule can be used: 

g 3 ^(a■n^^ an ^ , - rUt*''J 

3. yt/t — I 

and Var (z n ) for yU s l 

Higher moments, if they exist, may be found by a similar pro- 
cess utilizing the higher derivatives of f n ^-j_(s)„ 

A subject of great interest in any branching process is the 
probability of extinction of the species.,-' The Fundamental Theorem 
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of branching processes tells us that this probability is a fixed 
point of the generating function. The proof given follows that of 
Bharucha-Reid [2]. First, we note that if "a" is the probability of 
extinction, then aslim p^s lim f-fO). 

Theorem 3 . The probability of extinction, a, is the smallest pos- 
itive number, p, such that : 



Further, a = 1 if and only if JH £ J,, 

Proof : We see ^t hat a : 0 if and only if p 0 - 0 and f(l) = 1, so 

that 0-0-1 as it should be. 



Now f(s) is a power series in s with positive coefficients and hence 
strictly increasing for real arguments. If 0<a<b, then f(a)^f(b). 
Also, since p no >0 for all n, then a2 0. 

Let p20 be such that p= f(p). It follows that 




To show that a = f(a), we note that since p nQ = f u (0). 



Pn*,,o 




p jo rf| 0 )i f(p) = p 



Assuming p no ^ p, then f n (0)^ p and 



Pm 1,0 ' Pnt(fo) - Rp) s p 



So, by induction, p no £p for all n„ 



Then lim p no £p and a£p. 



It remains to show that a = 1 i'f and only if j t 



Assume first that a = 1 and show that ^u. ^ 1 
Now, it may be noted that : 

f 1 (s) 2 0, f"(s)Z0, so f(s) is convex in (0,1) 

f(o) - v Q yo , f(l) = 1 



f(s)y»s for s in (0,1) since 1 isthe minimum positive number 
satisfying s = f(s). 

Therefore, 1 - f(s)^,l~s and. 
l“f(s) 



1— s 



is bounded by 1 and is 
monotonically increasing 
with s for s in (0,1). 




Now assume Then, since f(s) is convex, f'(s) is either 
constant in (0,1) or strictly increasing with s„ By assumption, 
p-^ ^ 1 so P Q >0 and f'(s)^ 1 in either case for s in (0,1)„ 

By the mean value theorem, p IU“ f(s) - f wkei rC sO(<l 



Then, /-p(s) = p l (x)(l"S) and 

/- P(S) ^ 1-5 

f ts) > S ■ Por 04S*i 

Therefore, a = 1. 

Using theorem 3, it is a simple matter to determine the limits 
of the p nr for r ^ 0 o 

Theorem ii . A discrete branching process either dies out or becomes 

extremely large, i.e OJ lim p nr = 0 for r £ 0„ 

n->cP 



= 9 - 



Proof : Let t be a real number such that 04 t< 1 and recall the 



iteration property of f n (t)„ 

Case I. /A 4$ Then f(t)^t (see proof of theorem 3), and 
f(t)* f 2 (t)<f . . . . . . <1. 

Therefore, lim f r (t) = 1 s a, 

* oo n ' 

Case II 0/ u>J and t<£ a (p Q >0)» Then f(t)>t since f(s) is 




the limit of the last expression gives: 

oO 

lim f (t) r lim p + lira 3p„,„t r = a from the above three cases, 

n-»oo n *no ^ y ^ 

oo 

Therefore, ^Ildin p nr t r = 0° Since the series converges uniformly 

in an interval one may pass to the limit before summing, 
oo 

Therefore 3 lim p nT ,t r « 0 and lim p nT , = 0 for r 4 0 by the 

nr nr 

uniqueness of power series. 

Theorems 3 and h give us some interesting information. If we 
are working with an atomic pile, the probability is one that it will 
die out unless it is supercritical; then the probability that it 
will blow up is 1 - a. 
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Lotka [l7] found that the probability law of the number of sons 
a man might have was closely approximated by a type of geometric 
distribution, namely, p Q = 0.1*825, p k s (0.2126) (0„5893) k “ 1 ki- 1 
based on a 1920 census. * It follows that: 

f ( s) s 0.1*625 + 2^2126s f'(s)j * lo 2622 

1-0. 9893s 1 s-1 

Solving f(s) s s, it is found that a s 0.6226. Based on these values 
the probability of extinction of a family name is (0.6226) n where 
n is the number of male members of the family at the present time. 
This formula follows from the fact than, if there are r members of 
the zeroth generation then the generating function for the first 
generation is (f(s)) r . Then f n (s) * f n _^((f(s)) r ) j if f(a) sa, then 
(f(a)) r s a r and a r is the probability of extinction. 
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MEAN GREATER THAN ONE 



In this section we will consider the case wherey^ is strictly 
greater than one and the family has a non-zero probability of 



become extremely large, we will work with a normed random variable,, 
z n 

Let w n = — jf“ 0 The mean and variance of w n are easily deter- 
mine d 0 



and Bharucha-Reid [2^ „ We will prove first, the convergence in 
probability and then the convergence with probability one. 

Theorem If ju> 1, the random variables, w n , converge in probability 
to a random variable w„ 

Proof ,, Let n, m, be integers and n>m„ 



surviving indefinitely. Since we have determined that z n may 



= £(“-) = ~ & (Z n ) r 1 





Treatment of the convergence of w may be found in Harris f 10 ] 



£ Un|Z w »r) - X J [coefficient of s3 in(f n _ m (s) ) r ] 
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E(z n Um = r) = r[ F„. m (sjf"' f n '.„(s) 






5*1 



M- 






= r /X 

OO 

E(2 M 2 m )=| 6 J^ k E ( r2 n| 2 m a d 

- 2: p rV’^ 

» ju.* mm Ehk) 



E (w n w M ) ~ £ ! 2^, Zm) 






- £k J | 

f? [(UVr^mP] “ E )- E (u^i 



- _£L_ I J L 

" /SyU \ fjT ^ 



h-+<& 



™ Gju^-Wmf] - 0 



Hence the random variables converge in mean square and by a theorem 
of Kolmogorov <13, p.3U, ij , in probability 0 

Convergence with probability one may be proven with the aid of 
Doob»s fUj Martingale convergence theorem 0 



- 13 “ 



A discrete parameter martingale is a stochastic process such that? 



a ) E ( j y | ) 4 e>o for all n 

^ ^ ^n-tl I ^n» ° o ° s yi) ~ yn 

The theorem may be stated as follows,, Let (y n , n'ZO) be a discrete 

parameter martingale. Then E( |y 0 | ) £ E(|y-jJ ) £ E( y 2 )4 . . . 

If ^lim E( y n ! ) = k ^ oo f then lim y n s y* exixts with probability one 
and E(|y*i) £ k. 

Considering our random variable w n , E(w n ) 3 1 for all n. 




Hence the (w n ) form a discrete parameter martingale and we have 
Theorem Ua . If ju> 1, the random variables, w n , converge with probability 
one to a random variable w. 

Since w n converges to w, we may discuss the probability distri- 
bution of w. 

f 0 * 5 

Let G n (u) s P(w n £u), lp n (s) = E( e w n s ) ■ ^ e su dG n (u), 

pCiC 

G(u) s P(w^ u) , ip(s) 5 ^ e su dG(u ) 0 The distribution G(u) can be 
called the asymptotic distribution branching from f(s)„ There is a 
very interesting functional relation between the moment generating 
function (P(s) of w and the probability generating function f(s) as 
given by the following theorem. 
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Theorem 5 . (^) {^U ^) - f(Cf( S)) 



Proof : 



= fn(e£"j 



**♦.(«*) = FCf J^"l] 

- fCVnfs)} 



6 (s) , therefore, 

i 

Jr~CV^>) 

= C* su dGM = 4>(s) 

Upon taking the limit of (yUS)= f [ CP/) ($)") the desired result 

is obtained. 

It also follows that the k^h moment of w exists and is finite 
if and only if the k^h moment of exists and is finite,, 

Since G n (0) s P(w n £ 0) = p no , then lim G n (0) s a and P(w = 0)sa„ 
The following theorem will be used later. 

4 - 

Theorem 6 . Let G-j_(u), GgCu) be distributions having equal first 

# v. 

moments and finite second moments such that their characteristic 
functions it) and 1^2 (it) satisfy f [^, (^]J r = 1,2 

Then G-^(u) ~ G 2 (u) 



i <”i 6 n ( s) - 

h-><=o 

!iyr> LP„ls) = 

Kl-»eO n 
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Proof : (/? (it) r | + £ [ L i/+ oft") 

V--I V! ' 

for -t small by a special form of McLaurin’s theorem for small 
values of !tj . See [Cramer 3, p»27j . Then, 

i if 2 ut>i = 1 1 + f^tur + 0 (t 3 ) 

, tiil n * 

a. p) , , , . 

where /« ^(lt) = C**' 

| tf, U*) - tfa(tt)| = ^‘Jt^+ott 3 )! 

Therefore, J t/, (it) - Lfa ( t & ) | ~ f^3 (t) idksirg ^(i)< ' 

Then f tP, Ufyx) - (Jal/iU)} = | f Ctfj P C^(U)j| by theorem 5 

— yU j L^j (t*i * 1^)1 since when |$UJ 

Hence ? ^ (p<i ) j — ' t 

and !(S(t)| implying that JiA - j 

Butytt>i so ^?(t) must be identically zero 0 
Therefore j If, ( L t ) - lP A 1 1 * ■ ) ! ~0 

and CPj ( it) H (it) 

Finally G^(u) s G 2 (u) by the uniqueness of moment generating functions 0 
In order to determine some properties of G(u), it will be con- 
venient to work with a closely related function H(u) which is defined 

as the asymptotic distribution branching from k(s) where 
i^(s) - fCsQ-qj+ofl -a 



v 
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To demonstrate that k(s) is a generating function use the 



fact thatytol implies f(a) ~ a4 1. Then 

Y< o)s f'? 1 ' ' a = O 

I < f I ) = ^ — - j 

S — cx 



Clearly k(s) is monotonically increasing since,, if t^s, then 



KU) 



f : t ^ f [s(j-a)-tQ]-q _ | 

/-a ” ' ~ ~ K U J 



/-a. 



Therefore k(s) is a probability generating function. 

It is interesting to note the following properties. 



~ f ‘ EsQ-Qj+alO-a) 

/ * £k* 



f* [s(i-a)f a] 



/<’(/) - f'(!) ^ ju. 

f"[S{l-»)+a](|-a) 

K" III - (S A t/x^-yull l-a) 

K“(i) + K'(i) - CX'O )] 1 - s^d-al-al^-yu) 











r 



The coefficient of s^ is (l-a)^“^( p + P Ko. P G n + . , 

fk Tk*-» ' n / iKth 

and 9| = P, + + 3 * . , , = f'(Ci), 

Let l^(S) s l ' S J H bJ <> By a line of argument similar to 

the proof of theorem 5, it can be shown that ljp(y^5) = K C l S)J , 
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To demonstrate that k(s) is a generating function use the 



fact that JU>1 implies f(a) z a4 1. Then 

k' (o)= ----- — = O 

- I — CA 

Khl = h -' 1 ~ = j 

I 



Clearly k(s) is monotonically increasing since, 

^ ffslj-Ql-tal-a 

MU “ * /-a “ /-a 



if t/.s, then 

= K(s) 



Therefore k(s) is a probability generating function 0 

It is interesting to note the following properties. 

/{’(<) = £_Cs(Haj^] (i-q) _ f'‘[s(i'-ahd] 

! ^ (X 

K 1 ( < ) - r 1 ( I ) - yU. 

!< ’ 'S 1= f*Csd-ft|ta](l-a) 

«" (l) 5 ( S’ 1 +fJ?-JUL ) l I ' fi) 



K " ( 1 ) t K‘(I)- [Khir - 



CO 



Kls) = -/“'s-(z.P r !SH-a)+a) r 










The coefficient of s k is (l-a) K ~ 1 ( p i P K& + ■> « f tf u P O. f . » 

IK >K*-l » n f IKfh 

and q, = p ( f ap A a t3 + . . „ = 

Let \p (§) s 4 *" S (J Ljj , ( yj 0 By a line of argument similar to 

the proof of theorem 5, it can be shown that ^(^£^5) = K C ( S)3 , 
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We now proceed to determine a representation for the random 



variable w c First we determine the relationship between [fj ( s) and 

lfi( s). 



Theorem 7 e 



^(5) " 



/- a. 



Proof: Let W,(>)s. (PCti-a-iSl - O- 

/ 

Then (S) a MLCb&J5lLLl &). _ ... - 

/-CL 






and !///(«)= 



Therefore ! ~Pt and ^ have equal 



first moments,, Since (p!‘ 1$)' 



i = 0 






< <=^ 






IjA and l|/ have finite second morents c 

walls) » (PC(t-ft)Syu.i-fl. . f [ - a 

1 ' /-a /-a 

s £C{i~a) 4Ms)+&3 - a 

/-CL 



ty, (yUS) - KC Wl (S)] 



Therefore!^ and $ satisfy the same functional relationship with k(s)„ 
Therefore by theorem 6, 4>{S) = l|> f (s). 

Now we can determine the relationship between H(u) and G(u) c 
Theorem 8 . /^( U ) a, 6 

/ - a 

Proofs From theorem 7, IP (s) - (POi-QjSj-Q , 

/-a 



By definition of generating functions 



and 



/- a 



r°° uU-als 1 « . 

1 € d S ( u) “ ^ 

Jq 



y (si - b su dHM 

1 - u>C(i-a./s1-a 

J ' /-a 



Hence 
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^°°(i-a)e SU J Hl u l = \ e W, ‘ o1 - d <5(^-0. 

J 0 ” ° 

Let f(u) be defined as follows: £ ( u) - 0 f 0>t U £ 0 

£ (u) “ Cl tc* 'P'O 

Then ©■ - \ d £(o) as f 
Jq ** J 

Therefore a is the Laplace transform of £ (.0 ,). 

Hence j“f SW (| -<X)JW(W 3 6 (l^j - S^JcM 

and = f e sl 'd[<5f-;!''a)"£WJ 

v 0 

Therefore fj-fl) H (u) = 6 (* } ^a j " £M 

and HM « 7»o" [ ^(]^a )“ a ] for 

Clearly H(u) = 0 and H(«>o) - 1 as it should c 
To put the relationship of H(u) and G(u) in terms of random 
variables, we prove the following theorem G 

Theorem % The random variable w is distributed as the product of 

two indepei dent random variables, w Q and w', where P(w Q s o) = a, 

P(w„ s -i — ) s i_a and w* has the distribution H(u) 0 
u 1-a 

proof - e(£ hu/s J s a + (i-A) i-a J HM 

s a + li-ojC e 5 "J H lv(i-<4 

J O 

~ j o <* SJ J G ( u ) 

* ip Is) 
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is given with some 



So, P( u>(i-a) * U \w>o) - H(uJ. 

For later use a theorem due to Harris j lOj 
remarks. Let J = lo^ L_) * lo^-^ ) . If q ± ■ 0, let J^ob. 



Theorem 10 . If y<oo, Re(s )£o , s i 0, then IfMsJ = — - f Mo(s) 



where M(s) is continuous for s / 0, M(yWs) = M(s); M Q (s) -=■ 0 ? , Isl-* 00 

M(-s) = M(s). 

Remarks: (a) Under the conditions of the theorem M(s) is real and 
poaitive when s is real and negative. 

(b) If p(z, r j<io o, the r ^* 1 derivative of ^(s) satisfies 

l|As) - °(~rTrj 



(c) If f = oo , lp(s) and as many derivatives as exist approach 
zero exponentially as |s|-»<?o. 

We will now determine some of the properties of G(u) ; however, 
it will be convenient to prove theorems for H(u) and then interpret 
the results in terms of G(u). Let h(u) = h'(u). 

Theorem 11. H(u) is absolutely continuous. 

Proof ; Cramer [3] breaks up an arbitrary cumulative distribution 
function F(x) into three pieces: 

Fix) = Q| F, (x) 4 Q-, F A M 4 F 3 lx) 

where O ( f Ct^ f Q 3 » 1 0.^ Z 0 } L = l ; 3 

F-^x) is absolutely continuous. 

F 2 (x) is a step-function and is equal to the sum of the saltuses 
of F(x) at all discontinuities which are less than or 
equal to x. 

F^(x), the "singular" component, is a continuous function having, 

almost everywhere, a derivative equal to zero. 
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/ it x , 

O 2 l|X it) - /'■ p £* f is the sum of an absosutely convergent 
hVo * n 

trigonometric series and is thus an almost periodic function which 
comes as close to as we please for arbitrarily large values of t„ 
By theorem 10, lim If (t) - 0„ Therefore, for H(u), a 9 " 0 and H(u) 

£-»oo c. 



is continuous. 



Let h m (u) = i d{1 IPM Jt 

< rr\ 

dv - e dt u- (/ Hit ) 



Let 



v = -ft e 



‘it u 



h m (u) = 



Amu 






du - d [. W 

w ~ it { 






+ j 



i 

Am u 



vr 

I In/ V ,, , * 



A (WU) 

At 



Jt 



(a) h m W) = - I I/Me 



Examine each term of (A) as m oo 



1 j/l-ur)e 



I { 



"V 1 *" Jipm 

Am u jt 

rn ® 1 



First 



- Lmu 



MUm) 

Irni 



-trr.y 



T- + M 0 (s)j e‘ 

i Op-j.) 

iml j 



by theorem 10. 

- Lmu 



As we pass to the limit, this term is a function of continuous 
functions; hence, it must be continuous^ Similarly the second term 
leads to the expression [j> ( ~ t ) t — ; ml Milinj f ^H ay j] e 

which converges to a continuous function,, 

The last term In (A) Is C f — ^ Jt . 

qLttl u dt 

Since j C p/y) dx ■ ^ \ dtf , it follows that 

I J ! “ j 

I r f±“ im> j t j cC jt 

J ~ into Jt ~ J.no Ainu Jt 



- Jmu Jt l dL 
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by corollary b to theorem 10 
for some constant k 0 




'• =2niu 

- •Y) 




6 K 0 




Hence this term approaches zero as 0 Therefore, the continuous 

functions h m (u) converge uniformly in some interval, u-^ ^ u ^ Ug, 

to a continuous function h(u)„ 

M ~itu, 

H(u a ) ~ H(u,) = \',ry) f ty( L t) cjt 

J ^ - 3 m t 





e~ Ltu dudt 
t f(U) di 



~ Urn 



r u A 

Jy, ^Moll 



since h m (u) converges uniformly 



to h(u) the limit may be taken inside the integral and 



HM-HM -ztfhhlJu 

Since H(u) is continuous at zero and we have shown 'bsolute con- 
tinuity elsewhere, it is absolutely continuous on an interval in- 
cluding the point zero. 

In case£ Z^cO and $+ K-l , integration by parts of (A) and 

corollary b to theorem 10 shows that the first r derivatives of 

h(u) are continuous if u / 0„ • The integral expression for h(u) 

( r \ 

in terms of tp (it) shows that implies h v '(u) is continuous 

at zeroo This implies that G(u) " a Jh glvjdv for u>0j g(v) 
is continuous for v ^ 0. If then g( r ) (u) is continuous for u^O 

provided r<v^/<-j and is continuous for u s 0 provided r< y- J „ 
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Corollary to theorem 11 , p nr = P(z n s r)->0 uniformly in r, rZl 
as nsoo, 




Since G n (u)->G(u) uniformly for u>0, G(u) is uniformly continuous 
for u>0 and G(u) is right continuous at u s 0, the corollary is 
proven c 



It is interesting to note that Yaglom [22] treated this subject 
in a manner similar to I!arris 0 Yaglom defined 



°h s P to 

~ Q ft ~°h w r) 

H„(v) - P(y n “Y|y n >o) 



^(yJ “ limn H n { y) 

f e Liy JH( y) and proved the following theorem 0 

00 

Theorem 1$ 0 If yU>l , the variance of z ^ is not zero, and f"(l)4co 
then " y dt] ^ ()-&/ id 



, ' In this paper we followed the general line of Harris' argument 
and defined k(s) such that i^iyUS) ~ — t&l ~9r — ( 



This led us to the conclusion that H(u) was the distribution function 



Zyi 

for (l-a)w which is the limit of (1-p no )-ir 
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MEAN NOT GREATER THAN ONE 



We have examined the case for the mean greater than one and 
now we turn to the situation which corresponds .to the problem of 
subcritical or critical atomic piles 0 From theorem 3, the prob- 
ability of extinction is one and our fissioning process will 

% 

eventually die out. 

When the mean is exactly one, the limiting distribution of the 
conditional random variable, y n , (which was defined in the last 
section to be the number of members of the n^ generation, given 
that this number is not zero) is the same for any discrete branching 
process. 

Let a n = l-p no j = ijf- = P(z n = r ’ z n i 0) 

Qq 

F (s) = ]>” p* s r the generating function of the p'* r 

n r=o m 




if the limit exists 



Uniy) - P (% ^ y| Yn yo] 



Theorem 16 0 



hm 

n+oo 



If /UL = I t f S (l) ^0 , f ' then 

H n iy) - H (y) - / - e~ Y y^o 



-0 y 

(The proof may be found in Yaglom [22] . ) 

Now we will consider the limiting distribution of F n (s) in 
the casevxZl. 



-2U- 



Theorem 17 „ If ! and f'^cO 


then 


li m p* - p* 

¥) oc » Y if 


,22. v 
>. P*-./ 


linn F (S) s F(s) 

h-*o« **. 

and FCfl5))=^t F(s) +(h/u) 


|sU 


F(i)-I r ‘CD = K ' 


•J 



Proof : 

F (s)=f P(z„=r |j n *o) s r 

r- 1 



CK> 




r- 1 



^:^aioi s r 

P(Zn*D) 



by definition of conditional 
probability 



— ?{7n-r ) r 

fr; Pl2„iO) 



since the summation is over r 0 




The following theorem may be found in Konig ^lU’o If 3(s) is analytic 

in a domain, R, such that @(s Q ) = s Q and G\s q ) ! = a41 for some 

s in R, then, in this domain — converges uniformly to 

a 

a continuous function D(s) such that D’_0 ! S)J -Oil ($) 

'*■ 

Dls.1-0 D'i %)-\ 

- 25 - 



f(s) is analytic in |s'< 1. Since f(l) s 1 and f’(l)< 1, (yU<l), 
the s Q of Konig's theorem is one» We can conclude that £- n (s)~ i 

a” 

converges uniformly to a continuous function D(s) such that 
, (2) D(f (s) ) =/fD(s), D(l) = 0 , D * (1) = 1. 



' F n (s) ~ li 



frs(s )- | \ 

M 



n 



l i 



Fyi IQ ) 



Lim F fs) -P( S) = llir 






(3) F(s) - 1 + D(s) K“^ by Kordg's theorem, since convergence of 



i S uniform, the limit of exists „ 

jll 

From (3), then, 0 ($) - 
From (2), Off 1st) — y* D(s) 



= JU.KI Fui-a 

FC(is)j - l i by (3) 

K 

„ , 4 mKCBsH) 

“ I + /A- FI s) ~*yW. 



And, finally, Ftf(s)] = yU F (s) f I !-yuJ 
From (3), P(l) ~ | -f DO) K~' 

— ) since £) ? I ) = 0 by(2) 





D'(s) 



k' 



s- 



by (2) 
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One of the more interesting questions in the subject of branching 
processes is the determination of the number of generations to 
extinction when yUPl 0 

Let m = mm In: z h „ :0j 

oo 

ns 



6(5) -f R(M»n)e" s 

f)~0 


OO 

nrb 


b r. | - p 
J n r rH»,o 


b 3 = 


OO 

e, - J b„e' K 

t\*o 





''V / p « p \ p' 

l 1 ntB.O j ^ 

r-v-n 7 

Po 



Theorem 1 8, The probabilities b n satisfy the recursion formula^ 



- i 



n+! 



r (l~b n ) 



Proof : b 



oc 



r- *n + ii = & 



KiAjO ~ 1 P, ( P n+I 0 S°) 

- •> p,(i-b„y 

jk. J 



— . J - f(t~b 0 } 

Theoren 19„ The moment generating function of N is given by 

= I + (£“ 0 ^i(s) 



Proof s 



S 



oO 



> l Pn.,,0 - Pj e 



ns 



r^o 



GO 


(< 


y 


OO 

.> ( f 


r.n 






mi? 


=■ J 


4* 


oO > 

/ S(»i+i) 

.<1 e - e 


+ ®n + i 






/)rO 








oO 




- I 


f 


J (e - i)e (i 

ns 6 


~ Pn-H, 



no 









cO 



/ + b e 



5 



nsO 



I + (e s -l) 9 ,( 5 ) 
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So 0(s) is determined once 0-^s) is known. If we let z = e s , 
( s ) a power series whose coefficients are successive iterates 
of f*(b) = 1 - f(l-b) since. b n4l = f*(b n ) = fj +1 (b 0 ), where f*(0) = 0 
** (o) =/a< i. (see Fatou [83.) A function of this sort is meromor- 

So 



phic with poles at s 

y 



- n log jll , n s 1^2^ 

8,151 v 'T$ e ‘ 



o o o 



V-re" ' ->-e 

converges everywhere except at the poles, (see Lattes [29] ) 



oO 



x ( s ) = x iS i is determined by X(us) " f*(X(s) ) ,x"(l)s X*(l)s 1 

is 7 7 

y 0 is found by X(y 0 ) s b 0 s 1 - p 0 or by using inverse functions 
X*l(f*(s)) s yj. X“-*-(s) and determining X =, ^-(b 0 ) o 

If our process is of the type in which the "parent" does not 
die, then the total number of individuals is of interest. If fA 4 1, 
the probability that this number is finite is equal to one by 
theorem ija. 

o° __ 

Let 2 n ', then p(S<cOJs| f or ^ £ j 

‘ t oo r 

Let q r = P{ 5 ~ r) ; 

f “C 

Theorem 20 . Ifyui 1, then g(s) = s f(g(s))„ 

n 

Proof ; Let ~ 2 h 2»-Q £or some n with probability one 

1*0 

when s,‘ 1 by theorem Ua. Therefore, S n converges to some random 

variable S with probability one. The random variable, Sn, represents 

the total number of individuals in the generations from the zero^ to 

* 

the generation. Let $ r - * ^Yn where hVz \> ? since 2 b ~ S . 

cO ^ it 5 (I 

Let S K { S =2^S r ' , where 5} c . = P. and 5 - P ! y r - r J 

The event y^ - n can be divided into n mutually exclusive events 
with probabilities : 

P(zi “ i)P(n-i individuals in the following k-1 generations) 



for i= 1,2, . 
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So ©(s) is determined once £^(s) is known„ If we let z s e s , 
g-^(s) is a power series whose coefficients are successive iterates 
of f x (b) ” 1 - f (l=b) since b n4l r f*(b n ) = f*,^ (b Q ) , where f*(0) " 0 
f*'(0) -/*< lo (see Fatou L8]„) A function of this sort is meromor- 



phic with poles at s 2 - n logyu. 5 n § 1,2,, 



So 



/ c \ « Ml. + XaX 
,(S! hjue* ;y e s 



'v* 3 - 0 



converges everywhere except at the poles D (see Lattes [29] ) 
oO 

X(s) s x^s*- is determined by X(yx s) = f*(X(s)) ,x"(l )2 X*(l)s 1. 

Lil 

y 0 is found by X(y 0 ) s b Q s 1 <= p 0 or by using inverse functions 
X=l(f^-(s)) s s) and determining X“L(b o ) 0 

If our process is of the type in which the "parent" does not 
die, then the total number of individuals is of interest 0 If - 1, 
the probability that this number is finite is equal to one by 
theorem l*a 0 

00 _ V 

Let S = J 7^5 then p{S<cP)^| f or ^ £ I 

ir ' z ° vo r 

Let Cjj, = P(S»r) ; for 15 *'■ 

t oC 

Theorem 20 «, If j\*L 1, then g(s) s s f(g(s)) 0 

n 

^ 2 h 2r"0 for some n with probability one 



Proof : Let 5 ~ 



when ,(/ £ 1 by theorem Ua 0 Therefore, S n converges to some random 
variable S with probability one, The random variable, represents 
the total number of individuals in the generations from the zeroth to 



the nt* 1 generation.. Let 



Let 



n 






r - S ^Vru where V r s ' > 2 r , since 
, where ^ c= p and - P(y^-r). 



The event y^ ~ n can be divided into n mutually exclusive events 
with probabilities: 

P(Z]_ s i)P(n-i individuals in the following k~l generations) 



for is l,2,o 



= 28 - 



Starting with one individual at time zero, the probability of pro- 
ducing n individuals in the next k generations is the coefficient 
of s n in G^(s) 0 

Therefore, q, - . P(yi s i)P(n-i in the next k-1 generations) 
K a n S 

7 i i 

jr; . 

= > p^( coefficient of s n ' =1 in (Gj^iCs)) 1 ) 

Y, 

~ \ p^( coefficient of s n in (sG^^Cs) )^-) 

= coefficient of s n in ( Pi( s< \>i( s ) )^) 

s coefficient of s n in ( ]> P^isG^^is))^) 



~ coefficient of s n in f (sGj^^s) ) 

Therefore, G k (s) = . q k n s n » f(sG k=1 (s)) 0 Clearly, g n (s) s sG n (s) 



ns o 



since P(y n = n) = P(Sn = nt 1) and g n (s) ® sf (g n ^]_(s) ) „ Taking the 
limit as n approaches infinity, we see that g(s) - sf(g(s)) 0 

Using theorem 20, the expected total number of members of a 



family is: :?'(?)- 4- ' qU'f] Z r cf/sH 9 ! ( Si ! 

■y r\ - ^ s ? r * * * 



Sr. J 



= i 4 •; <*'o) 

0 



Then, 



q " { i 






SlBilarly, p"fl) » Mj zM 



and the variance of S is; 



I 



= fU) - r/l :)-■ 



'* ' / 






{‘-mi 3 



y:/.. 






( -/A 4 (i- 
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f (i ) ( i-yu)(/U~ l) f 

( 1 -/^ 



5 '-//P 

S 3 - yU -^/U a 
( 1-yUj 3 
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ESTIMATION 



In the introduction we mentioned that the parameters, p r , must 
be empiraca .ly determined for the given amount and arrangement of 
atomic material. The method of maximum likelihood can be applied; 
however, it seems necessaiy to assume that certain random variables 
are observable. 

Let be the number of individuals in the m^ generation 
who have exactly r descendents in the (m4 l) s ^ generation. Then 

m 



oo 

2 rn = — 

f-0 



oo 



rz 



m-i, r 



Let S n 2, 



r-o r ~ 0 

Theorem 21 , The maximum likelihood estimates of p r and JtA based on 
observed values of for m'in (n41 observed generations) are 

p r - y a - 

Proof : In order to construct the likelihood function, first we 



determine the conditional distribution of z^,, r ■z 0,1, 



© © © 



given 






From our initial assumptions, each member of the generation 



has the same probability distribution associated with his progeny. 
Therefore the probability of k Q members having no descendents is 

P 0 ° | J<*| where [ ^ j xs the binomial coefficient and is equal to 



2 

_ -fr ■ „ Of the remaining z m - k D members, say z m - k, have 



one descendent. This probability must be 






w • 



Continuing 



in this manner we get the multinomial functions: 



^( z mo> z ml> 



o » « 



u mi> 



9 © O 



| z m-l ,0> Vl,l » • • •> z m-l,i> •••) 

co _ 



7 J TT 



7T fZr»r)j 

Y-0 
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. . . 1 


Z n-1,0> z n-l,l> • • *)P(z 


n-2.1 » 




II 


OO - 

-7 I 7 T P / ™~ lr 

yys J r r 


fr |Z„„)/ 


m-o 



Ir-to 

Due to the Markov properties of the branching process this is the 
joint distribution of the , m = 0,1,. „ . , r= 0,1, . . . 

The only factor of the likelihood function which depends on the p r is 



oo 



7 T p r 






oo \ 

the logarithm of this expression is | I ofy ( ) / 






To maximize this expression subject to the restraint, p — | 



the method of Lagrangian multipliers may be used. 

L=£ it p K -Jtp ) 

f s o m - o 
Then Z. — 



K=® 



2>P, 



r ” ^ ( 2 ™r(pH -/)) 

Y yw=o 



Seeing -|k - 0 , we get 2 Wh for r = °> 1 ^ 






But p s 0 5 summing the above expression over r gives: 

r - 0 oo n 

2 ]> 

Y-o ri\~o 

The double sum of the z^ is merely S n so ,)| = S n and 

✓v ■*— 

Pr = 






oO 



To find p. recall that yU - ^ > | rp^ and that IT - Z w+1 



r=o 



r=p 
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Using these facts^ it follows that: 



oO 



A 



Y ^ Zrr\r 

-Es£ 



S„ 



<\ oO 



rz { 



mr 



— rn-o r= 0 



— 



s„ 






s„ 



P = 



— Sn+j I 
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ION FLUCTUATION 



The number of ion pairs produced by a fast primary particle 
passing through an absorbing medium is a branching stochastic process 
(Moyal [18] ). In this process the primary particle, which may be 
a gamma ray, has a certain charge, mass, and velocity or a given 
frequency; the thickness and atomic properties of the absorbing 
medium are known. The primary interacts with the atoms in the 
absorber, freeing electrons by the Compton effect. These electrons 
are the first generation and each may ionize another atom or atoms 
depending on their energy; secondary electrons may ionize more 
atoms, etc. The process ends when all the freed electrons are too 
slow to produce further ionization. This results in some number of 
free electrons and, thereby, an equal number of ion pairs. 

In determining the distribution of this number we shall make 
the following assumptions: 

1. Successive ionizing collisions are statistically independent. 

2. The total energy loss of the primary is very much smaller 
than its initial energy and, hence, the loss of primary energy may 
be neglected. Therefore, we are not treating slow primaries or 
very thick absorbers. We also assume that energy losses due to 
radiation, nuclear interactions, etc. are negligible. 

3. The absorbing medium is homogeneous. 

U. The recording of data is delayed sufficiently to allow 
termination of the process. Photoelectric effects in the gas or 
chamber walls is neglected. 

Unlike the branching processes treated in the previous sections 



the probability distribution of the number of progeny per individual 

in a given generation is dependent upon the generation,, Thus one 

may expect different methods to be used 0 Moyal has applied the 

Laplace transform method to solve the Kolmogorov equations of a 

Markov process and this is outlined below. 

Let t be the thickness of the absorber 

N = the number of absorber atoms per unit volume 

6"(E) = total cross section 

q = primary ionizing collision rate 

q n - P(n ion pairs are produced in a given collision) 

p_(t) = P( n ion pairs produced in a thickness t) 
n oo * 

M(s,t) = p n (t)e" ns , the generating function of p n (t) 

h=o 



Then (n 6(E) (dE^ dt^ is the probability of an energy loss between 
E and E + dE in an absorber of thickness dt„ The primary ionization 
rate is c] = N $ 0*te)d£ 

oO 

The - J ' and q(dt) - P( an ionizing collision in thickness dt) . 

n-o 

Since at least one ion pair is produced in a given collision, q 0 = 0» 
Based upon the above assumptions vre have 



(A > Pn(t| +*») = It Prt-jU*) 

JsO 



p n ( gt) =(I-C|($fc) r no + 1 o(H) where £„ 0 is the 

Kronecker delta and 0(£t) goes to zero faster than £t as £t-*0. 
By assumption 1 and equation (A) , 
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&M(s.t ) - 



mi M(S, t * it) - M(s, t) 
K gt-»o %t 



= [£0 Vr(M< s ^i)-l)M|S,tj 



Inrn -i. 



I - q(St) + ""HImM 

WsJ 






L h-i 



©O 



^ =nis.)hS,,[ e --,) 

Wi * n 



n*i 



Comparing this equation with "^jr~ ~ cy 5 it c l ear that 

M(s <t )= e xphtf( e -«-|)^]^ e « Rts) 

<0 = i J 

where <5 = c,t 



R(s) = 2 (^- 1 )^ 



The inversion of the Laplace transform, M(s,t), is given by: 



P n te> = 



1 

air L 



' c * nL qR(u) nu 
\ £ S 

t-rr£ 




(see Widder tap 



Evaluating this integral by the saddlepoint method (see Korn [l£]), 
using the first term only, gives: 

P„(Q) = -i-[an Q R"(s„)]" te e*p[q( RW- S„ R 1 (s„)] 

where S n is related to n by 

n=-QR'(s„) = <3 j> (Ke'* S ’7( c ii < ) 
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since 



R(sj 




R“(S„J = f 

K-i 



-S„K 



and c is a normalizing constant,, 

Slnce X P„!Q) - ! 

n-D 

QO t '/A | . 

c = Z[^' T QR"(S h )] exp [Q(R(S„)-S^R (s„ 

050 
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